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ABSTRACT 

We report on new radial velocity measurements of massive stars that are 
either suspected binaries or lacking prior observations. This is part of a survey 
to identify and characterize spectroscopic binaries among O-type stars with the 
goal of comparing the binary fraction of field and runaway stars with those in 
clusters and associations. We present orbits for HDE 308813, HD 152147, HD 
164536, BD-16°4826 and HDE 229232, Galactic O-type stars exhibiting single- 
lined spectroscopic variation. By fitting model spectra to our observed spectra we 
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obtain estimates for effective temperature, surface gravity, and rotational velocity. 
We compute orbital periods and velocity semiamplitudes for each system and note 
the lack of photometric variation for any system. These binaries probably appear 
single-lined because the companions are faint and because their orbital Doppler 
shifts are small compared to the width of the rotationally broadened lines of the 
primary. 

Subject headings: binaries: spectroscopic - stars: early-type - stars: fundamental 
parameters 



1. Introduction 



The multiplicity frequency of O-type stars in clusters or associations is large. ISana et al. 

(120111 ) find that ~ 44% of O-stars in open cluste rs are members of relatively short period 



spectroscopic binairies. iMason et al.l (119981 120091 ) considered the total numbers of spectro- 
scopic (short period) and angularly resolved binaries (long period) and derived an estimate 
of ~ 75% for O-stars with some companion, but this estimate must be a lower limit because 
the intermediate period range (years to decades) remains largely unexplored. Spectroscopic 
detection of companions to O stars presents many problems. The high luminosities of early- 
type stars can easily outshine spectroscopic companions with a difference of just a few tem- 
perature subtypes, and/or a single luminosity class. Some O stars are rapid rotators with 
broad lines that can hide the spectral features of a companion, even with high resolution 
spectra. Systems may also be seen at an unfavorable inclination to our line of sight, reducing 
the shifts of spectral features that would be observed. Some combination of these three sce- 
narios results in the detection of such systems as low-amplitude, single-lined spectroscopic 
binaries (SB1), or no detection of velocity variability at all. 

We present results on five such systems here that were observed as part of a larger group 
of O-type stars with little prior observati onal work and that a ppear with a binary designation 
of "U" (unknown) or "?" (uncertain) in IMason et al. Jl998h . The current status of each of 
these systems is either unknown spectroscopic variability (HDE 229232 and BD— 16°4826) or 
known variability but with insufficient data to compute an orbit (HDE 308813, HD152147, 
and HD 164536). With new data presented here (§2), we measure radial velocities, determine 
orbital elements, and compute stellar parameters for each system (§3), close ly following the 



proce dure used to analyze velocity non-variables in Paper I of this series (j Williams et al. 



20111 ). We discuss these first single-lined spectroscopic orbits and physical parameters for 
each system (§4). We conclude the paper with a discussion about the nature of the unseen 
companions (§5). 
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2. Observations 

Observations of HDE 308813, HD 152147, BD-16°4826, and HD 164536 were obtained 
from 2003 March 16 to 20 UT, 2004 January 13 to 29 UT, and May 27 to June 8 UT with 
the Cerro Tololo Inter- American Observatory (CTIO) 1.5 m telescope and RC spectrograph. 
We used grating $47 (831 lines per mm, 8000 A blaze wavelength in the first order Littrow 
configuration) in second order together with a BG39 or CUSO4 order blocking filter with a slit 
width of 2". The detector was a Loral 1200 x 800 CCD. This arrangement produced spectra 
covering the range 4058 - 4732 A with a FWHM resolving power of R = A/ A A ~ 2750 
and a wavelength calibration accuracy of ~4 km s _1 based on the rms scatter of fits to 
comparison spectra. The instrumental broadening from this setup, found from the FWHM 
of the comparison lines, was computed to be equivalent to a Doppler broadening of 109 
km s^ 1 . Exposure times were several minutes in order to reach a signal-to-noise ratio (S/N) 
exceeding 100 per pixel. Each observation was bracketed with HeAr comparison spectra for 
wavelength calibration, and numerous bias and flat field spectra were obtained each night. 

HDE 229232 was observed with the Kitt Peak National Observatory (KPNO) 2.1 m 
telescope during an observing run from 2008 November 15 to 21 UT. These observations 
made use of the Goldcam spectrograph with grating G47 (831 grooves mm -1 ) in second 
order with a CUSO4 order sorting filter. The detector was the T3KC CCD (a 3072 x 1024 
pixel array with 15 x 15 /im pixels), and the slit width was set to 1'.'3 resulting in spectra 
with a FWHM resolving power of R = A/AA = 3030. The exposure time for HDE 229232 
was 600 s, enabling a S/N of about 200 per pixel in the continuum. The wavelength range 
is 3942 - 5032 A with a wavelength calibration accuracy of ~5 km s" 1 based on the rms 
scatter of fits to comparison lines. Instrumental broadening for these data, again estimated 
from the FWHM of comparison lines, was found to be 99 km s _1 . 

All spectra were reduced using standard routines in IRAF§ and then rectified to unit 
continuum via fitting line-free regions. Next, all spectra for a specific target were transformed 
to a common heliocentric wavelength grid in log A increments and were placed into a two 
dimensional array as a function of spectrum wavelength and time. The subsequent velocity 
measurements and analysis of spectra may then be done on the entire array of spectra, rather 
than on each individual spectrum. 



3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Radial Velocities and Stellar Parameters 



To measure both velocities and stellar parameters, we began with an approximate model 
template for the star. A simple "by eye" fit to several prominent a bsorption features was 
made by interpolating in the OST AR2002 grid (ILanz fc Hubenvil2003f) for stars with T e g > 30 
kK and in the BSTAR2006 grid jLanz fc Hubenv!l2007h for stars cooler than 30 kK. 



Velocities for each spectrum for each object were then measured using the standard 
method of cross- cor relation , and the uncertainty in the resulting velocity was estimated 
using the method of IZuckerl (120031 ). The model was cross-correlated with each spectrum in 
the particular object's spectrum stack. Prior to cross-correlation, object dependent emission 
features and interstellar features were removed, including the broad diffuse interstellar band 
(DIB) near A4428 A. The radial velocities are listed in Tabled] with targets ordered by right 
ascension. The table lists the object name, Heliocentric Julian date of observation, orbital 
phase (see below), radial velocity, velocity uncertainty, and the observed minus calculated 
(O — C) value for each measurement. Note that the quoted velocity uncertainties do not 
include those uncertainties related to the wavelength calibration of the spectra. 

For each object, we ran fits of the radial velocity measurements f or eccentric orbits and 
for cir cular orbits using the nonlinear, least-squares fitting routine from lMorbey &: Brosterhus 
( 119741 ) . We note that we did not use any other period search algorithm because our data 
sets covered all of or significant fractions of each orbital period, so there was no ambiguity 
in the selection o f an in itial period estimate for the general orbital fitting scheme. Following 
Lucy fc Sweeney! ( Il97ll ). we computed a statistic based on the rms of the fit for the eccentric 
and circular cases, p = (rms e /rms c )^, where rms e is the rms for the eccentric orbital fit, 
rms c that for the circular orbital fit, and = (N — M), where N represents the number 
of observations, and M the number of parameters fit (6 for the eccentric case). Values of 
p > 0.05 (the 5% significance level) indicate that the eccentric orbit does not statistically im- 
prove the fit. For each of the orbits computed here, the eccentric fits provided no significant 
improvement in the residuals, so we adopted circular orbits. The circular orbital parameters 
of period, P, the time of maximum velocity, T , the systemic velocity for the system, 7, 
and the velocity semiamplitude, K\, are listed in Table [2j Also listed are the mass function, 
f(m), the projected orbital semi-major axis, aisinz, the rms of the fit, and the number of 
observations. The radial velocities and orbital solutions are plotted in Figures [T] through El 
as are the uncertainties associated with the wavelength calibration, <j\. Using the final orbit 
fits we employed a shift-and-add algorithm to the spectrum stack in order to create a higher 
S/N master spectrum. These spectra are plotted in Figure |6j 

Stellar parameters were estim ated by fitting the mas ter spectrum for e ach object with 
model templates from OSTAR2002 JLanz fc Hubenvll2003[ ) and BSTAR2006 flLanz fc Hubenv 
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20071 ). Initial values for v sin i were estimated by convolving the model template with the 
instrumental broadening (see §2) and with a projected rotational velocity (v sin i) function 
to obtain a broadening function making a rough match of one prominent non-hydrogen ab- 
sorption line. The choice of the absorption line used differ ed based on the temperat ure of the 
star. We adopted linear limb darkening coefficients from lWade fc Rucinskil (119851 ) . A least 
squares grid search routine was then employed to refine the estimates of T e g and log g from 
fits to several spectral features (listed for each individual object in the following section). 
These updated estimates for T eS and log g were then used to refine the v sin i value by 
comparison of the FWHM measurements of lines in the master spectrum with those from 
the model template broadened for a range of v sin i. The adopted T eS and v sin i values 
are means of individual line measurements with uncertainties computed from the standard 
deviation of the measurements of different lines. For log g, we used a weighted average where 
Balmer line wing fits are given twice as much weight in the calculation of means because of 
their sensitivity to pressure broadening. It should be noted that we cannot measure v sin % 
values less than roughly half our i nstrumental broadenin g, or about 50 km s _1 . Another lim- 



itation is that both OSTAR2002 flLanz & Hubenvil2003f ) and BSTAR2006 (ILanz & Hubenv 
20071 ) are based on plane-parallel approximations and solar helium abundances. These ap- 
proximations may be less accurate for very high luminosity stars, such as HD 152147, which 
generally have extended atmospheres and significant winds. Lower values for log g indicate 
the star is evolved, and such stars may also have enhanced helium abundances due to mixing 
and mass loss. The entire process of fitting stellar parameters was repeated once to ensure 
accuracy of the final results. 

The derived T e s, log g, and v sin i values are listed in Table [3j along with a spectral 
classification estimate made b y com paring our computed values for T e ff and log q w i th Ta bles 
4, 5, and 6 of iMartins et al.l (120051 ) for O- stars, and Table 2 from iBohm-Vitensd (119811 ) for 
the two targets on the border between O- and B-type stars, HDE 308813 and HD 152147. 
Column 2 in the table lists the spectral classifications for the primary component found 
in the literature. The final three columns in Table [3] list the parameters used to obtain a 
spectroscopic parallax distance estimate that is described in detail for each object in the 
next section. 



4. Discussion of Individual Objects 
4.1. HDE 308813 



The vari able velocity of HDE 30 8 813 a nd its membership in cluster IC 2944 were first 
suggested by iThackeray fc Wesselinki (119651 ) based on three observations spaced two years 
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apart. They suggested that HDE 308813 may be an SB 2 due to differe n ces in velocities from 
He II compared to other lines seen on their last plate. iHuang fc Gied ( 120061 ) also measured 
three velocities, over a tim e span of 4 d ays, a nd suggested the system was a possible single- 
lined spectroscopic binary. ISana et al.l ( 1201 ll ) gathered 11 velocities over time spans of days 
to years, and they were therefore unable to discern an orbital period or ephemeris fo r this sys- 
tem, although they did confirm the SB1 nature suggested by lHuang fc Gied (120061 ) . The ob- 
servations reported here have the benefit of being made over three separate epochs, the second 
one a 16 day span. Our computed orbital period of 6.340 d (T abled is therefore covere d twice 
in that time span. The rad ial velocity measurements fr om iThackeray &: Wesselinkl ( 1965) 



(-12 18, and -40 km s" 1 ). lArdeberg fc Maurice! (119771 ) (-10 km s _1 ), and IHuang fc Gies 
( 12006) (21.9, — 13.4, a nd — 17.6 km s _1 ) all span the range of velocities we observe. The data 



from IHuang fc Gied (120061 ) are overplotted in Figure [fl but were not used in computation 
of the orbit. Including any of the archival data increased the rms of the fit by at least a 
factor of two (presumably due to differences in the way the velocity was measured by these 
authors), so we only fit data reported in this work. 



Prior v sin % measurements are 186±11 km s 1 by iDaflon et al.l (120071 ) (using He I 
AA4026, 4387, 4471), 196±9 km s" 1 bvlHuang & Giesl (bood ) (who used the same lines plus 
Mg II A4481), and 197 km s _1 by iPenny et al.l (120041 ) (from UV lines in spectra obtained 
by the Hubble Space Telescope and Space Telescope Imaging Spectrograph). All of these 
measurements agree within la of our computed value of 204±10 km s^ 1 obtained by measur- 
ing He I A4143, 4387, 4471. T eS and log g were estimated using H5, H7, He I AA4143, 4387 
and H e II A4686. Based on the values computed above and the calibration of iBohm- Vitense 
( 1l98ll ). t he spectral cla ssification of HDE 308813 is B0 V. This is very close to the classifi- 
cation bv lSchildl fjl970h of 09.5 V. 



W e can estimate the distance to HDE 308813 via spectroscopic parallax. IThackeray fc Wesselink 
f)l965h measured V = 9.28 arid (B- V) = 0.04 for HDE 308813. We adopt (B -V) n = -0.26 
for a B0 V star from IWegnerl (119941 ) and My = —3.90 for an 09.5 V star from iMartins et al. 
(120051 ). Assuming Ry = 3.1 leads to a distance of 2.82 kpc, and this value is listed in Column 
9 of Table El This distance puts HDE 308813 a bit further than the distance for IC 2944 of 
2.2 kpc from iTovmassian et al.l (119981 ) . but it may be a member of a more distant association 
(see their Table 5). It should also be noted that the difference in the distances ma y be due to 
the un certainties in our estimation of (B — V), (B — V)o, My, and Ry. However, iDias et al. 
(120021 ) report an average cluster velocity for IC 2944 of — 4.8±5.1 km s _1 . The systemic 
velocity from our orbital fit for HDE 308813 is 8.5 km s" 1 , 2.6<r from the cluster average, 
making cluster membership unlikely. 
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4.2. HD 152147 



HD 152147 lies in the vicinity of the Galactic open cluster NGC 6231, but it is not 
l ikely to be a m ember. A number of velocities have been previously published, five by IStruve 
(119441 ). one by IConti et al.l (119771 ). and 10 by iLevato et al.l (119881 ). who were also the first 
to claim velocity variability. We fit an orbit to all these literature velocities and our own 15 
measurements with equal weight and obtain the parameters listed in Table [2] with the orbit 
shown in Figure [2j 



Conti fc Ebbetsl (119771 ) measured v sin i = 80 km s 1 by comparing H/y, He I AA4387, 4471, 



and He II A4541 with calibrated standard stars. We measure a rotational velocity of 150±28 
km s _1 from fits to Si IV A4088, He I AA4143, 4471, 4387, and He II 4686. Their spectra were 
obtained at higher resolution than our spectra, so our measurement should be taken as an 
upper limit. The estimates for T e g and log g were made from fits to H7, He I AA4143, 4387, 



and H e II A4686. Using these numbers in concert with the calibration in iMartins et al. 



(120051 ). we estimate the spectr al classification for HD 152147 to be 09.5 I, very close to the 
spectral classification given by IWalbornl (119721 ). 09.7 lb. 



To estim ate the spec t roscop ic parallax, we began with V = 7. 2 3 and (B — V) — 0.37 



measured by ISchild et al.l (119691 ) from seven observations. IWegnerl (119941 ) gives a value of 
(B — V)o = —0.24 fo r an 09.5 I star. Fur ther assuming a standard value of Ry = 3.1 and 
using My = —6.28 in IMartins et al.l (120051 ) for an 09.5 I star, leads to a calculated distance 
to HD 152147 of 2.11 kpc. This value does not agre e well with the recent distance estimate 
for HD 152147 of 1.64±0.24 kpc from iMegier et al.l (120091 ) that is based on interstellar line 
strength. Nor do es this value match with the distance estimate to NGC 6231 given by 
Sana et al.l ( 120061 ) of 1.64 kpc. This may mean our estimate for Ry is inaccurate for the 
line of sight to HD 152147, but this also may be evide n ce su ggesting; that HD 152147 is 
not a member of NGC 6231. For example, iHumphreysl (119781 ) places this star in the Sco 
OBI association a t a dis tance of 1.91 kpc, a distance closer to our estimate. Furthermore, 
Kharchenko et al.l ( 120051 ) reports that NGC 6231 has an average cluster radial velocity of 
-27.28±2.98 km s _1 , far different than our systemic velocity from the orbital fit of — 40.1±0.7 
km s -1 . 



4.3. HD 164536 



The first radial velocity measurement published is by iHayfordl (119321 ) who reports a 
value of 6 km s -1 from five spectral lines in a single spectrum. Other radial velocities are 
averages from several, typically unspecified, nights: —11.1 ± 3.2 km s _1 from five plates 
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report ed by iNeubauerl (119431 ) and —10 km s 1 from three plates reported by iFeast et al. 
(119551 ) . None of these velocities were used in the derivation of the orbital parameters listed 
in Table H 



HP 164536 is one of th e brightest members of the cluster NGC 6530 ((Walker 1957; 

Boggs fc Bohm-Vitensdll989l ) . A proper motion study was made of NGC 6530 by Ivan Altena fc Jones 
( 19721 ). who claim that HD 164536 is not a member. Howev er, its cluster membership 



was reestablished in a series of pa pe rs based on IUE spe ctra (IBohm-Vitense et al.l Il984 



Boggs fc Bohm-Vitensd Il989l . Il990l ) . lAfricano et al.l ( 119751 ) claimed detection of a nearby 
companion during a lunar occultation observation at p = 0'.'174 ± 0'.'003 arcsec, 9 = 203°, 
and Am = 2.2 ±0.4 mag. Subsequent follow up observations that had the ability to see such 
a companion ( Masonlll996l ; iTurner et al.ll2008l ) did not resolve this close companion. This 
star is the brig htest of a trapezium (4 star) system liste d in the Washington Double Star 
Catalog (WDS; iMason et al.ll200ll ). ITurner et al.l (120081 ) resolve the closest component at 
l'/7 distant, with other members outside of their field of view. This close c omponent is too 
faint (V = 12.4 according to the WDS) compared to HD 164536 (V = 7.1; iHeig et al.llioooh 
to influence our spectra. 

To derive the effective temperature and log g, B.5, H7, He I A4471, and He II AA4541, 4686 
were used. The projected rot a tiona l velocity was computed from observations of the He II 
lines. iBoggs fc Bohm-Vitense Jl989l ) arrive at an effective temperature of 32 kK determined 
from optical photometry and spectral type of 09 V from their analysis of IUE spectra. The 
value derived here of 37.4 ± 0.9 kK is > 3a higher, i ndicat ing an 06.5 V estimate of the 
spectral classification from Table 4 in iMartins et al.l (120051). This classi fi cation is also at 
odds with the optical classification reported by iMacConnell fc Bidelmanl ( 119761 ) of 09 III, 
although the spe ctrum is noted as being overl apped with another star, making the classifi- 
cation uncertain. iHouk fc Smith-Moord (119881 ) find a spectral typ e of 07/08, closer to our 
estim ate. G iven the uncertainties in the s pectral classifications in IMacConnell fc Bidelman 
( 119761 ) and IBoggs &: Bohm-Vitensd (119 89), we think our estimate of 06.5 V matches well 
with that of IHouk &: Smith-Moore ( 1988 ). and is currently the best estimate. 



We estimate the distance to HD 164536 by taking the Tycho B and V magnitudes 
dHog et all lioooh transformed to Johns on magnitude s, an My = —4.77 for an 06.5 V 



( IMartins et al. 



20051 ). (B - V) = -0.29 ( IWegnerlll994h . and ass ume R v = 3.1. This results 



in a distance of 1.78 kpc, equal to the 1.78 k pc estimate given by Ivan Altena fc Jonesl (Il972[ ) 
and larger than the 1.33 kp c distanc e from iKharchenko et al.l ( 120051 ). Using Ry = 3.45 for 



HD 164536, as reported in IWegnerl (120031 ) . the distance drops to 1.71 kpc. Both of the 
distances computed in this work agree within uncertainties with the distance reported by 
van den Ancker et al.l ( 119971 ) of 1.8±0.2 kpc. The systemic velocity component of HD 164536 
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found here is — 4.9=1=1 .1 km s 1 ; this value is 1.3a from the cluster radial velocity of NGC 
6530 listed by lBarbier-Brossat fc Figonl (120001 ) of — 13.32±6.34 km s _1 . This appears to offer 



tentative support for the inclusion of HD 164536 in NGC 6530. 



4.4. BD-16 4826 



No earlier velocity measurements were found in the literature for BD — 16°4826, making 
our 11 velocities and SB1 fit the first published. However, additional observations are desir- 
able because our observing span of 12 d is somewhat smaller than the 16 d orbital period. 
Consequently, the orbital parameters should be considered as preliminary values. 

BD— 16°4826 is a member of the cluster NGC 6618 which is embedded in M 17. Esti- 
mates for effective temperature and log g were made from fits to B.5, H7, He II AA4200, 4541, 4686 
and He I AA4387,4471. We used He II AA4541,4686 to derive the v sin i value listed in Ta- 
ble El There was very large scatter in our attempts to fit the spectrum of BD— 16°4826, 
as evi denced by the highest unc ertainties for stellar parameters in this work. Listed as an 
05 in iHiltner fc Johnson! (11956), the spectral classifica tion of this object has varied from 
03 V based on optical spectra (IGarmany &: Vaccalll99ll ) to any of four spectral types rang- 
ing from 05 V to 09 III that were based upon IR data by iPovich et al.l (120090 . The best 
estimators we have available as luminosity criteria are the wings of the Balmer lines. These 
definitively rule out the possibility of BD— 16°4826 being evolved (see Fig. [6] and Table |3|), 
and we estimate the spectral classification for BD— 16°4826 to be 05.5 V. 



The most recent optical photometry is reported by lOgura fc Ishidal (119761). They find, 
from two observations, V = 9.85 and (B — V) = 0.94. Using My = —5.07 (IMartins et al. 
20051 ). (B — V)o = —0.30 (jWegnerl 1 19941 ). and a standard value Ry = 3.1, we arrive at a 



spectroscopic distance to BD— 16°4826 of 1.60 kpc. The distance to N GC 6618 and M 17 



is a topic of some debate. The recent IR study by IPovich et al.l ( 20091) lists two distances 



1.6 ± 0.3 kpc from iNielbock et al.l (boOlh and 2.1 ± . 2 kpc from 



The value derived here agrees with the INielbock et al.l (120011 ) estimate. Our c omputed value 



Hoffmeister et al.l (120081 ). 



i s also within the uncertainty of the lower limit on the distance reported by lHanson et al. 
(119971 ) of 1.3i°;2 kpc based on extinction corrected IR magnitudes of O-stars. The systemic 



velocity found here for BD— 1 6°4826 of 11.0 ± 1.0 km s 1 is 2.2a different than the cluster 
velocity of — 17 ± 13 km s _1 (IKharchenko et al.ll2005l ). This discrepancy in velocities hints 
that it is less likely that BD-16°4826 is actually a member of NGC 6618. 
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4.5. HDE 229232 

There is very little in the literature about HDE 229232. Like with BD-16°4826, our 
seven velocities are the first published. However, with only seven data points, our orbit 
should be considered preliminary. For example, it is possible that the system has a signif- 
icant eccentricity and longer period, so that our data record the variations around a single 
periastron point. Our estimate that the star is an 05 V with a variable velocity indicates 
that this is a massive binary that deserves additional observational work. 

To estimate effective temperature and log g, we used the B.S wings, He I A4026, and 
He II AA4200,4541. Only the helium lines were used to com pute the y sin i . Our estimated 
spectral classification of 05 V is very close to that listed by I Walbornl (119731 ) of 04 Vn((f)). 



Using V = 9.525 and (B-V) = 0.82 fro mlHiltner fc Johnson! (119561 ). My = -5.21 from 
Martins et aD tood ) and (B - V) Q = -0.30 ( IWegner!ll994l ) (both values for an 05 V star), 
and the standard Ry = 3.1, the distance to HDE 229232 is 1.79 kpc. The galactic coordinates 
of HDE 229232 are I = 77?40 and b = 0?93, looking down the Cygnus arm, where many 
associations wit h a range in distances from 0.75 to 2.32 kpc appear more or less in the same 



part of the sky (jUyamker et al. 



2001 



, and re ferences therein). It lies in the direction of the 



Cyg OB8 association that iHumphreysl (119781 ) places at a distance of 2.29 kpc. The systemic 
velocity from the orbital fit of HDE 229232 is -46.6±0.8 k m s" 1 . SeFNik fc Damb3 ( 20091) 
report a velocity for Cyg OB8 from velocities in the work of iBarbier-Brossat Sz Figo nl (hoooh 
of — 21±11 km s -1 . This puts our estimate of the systemic velocity 2.3a from the mean value 
for Cyg OB8, perhaps lending support that HDE 229232 is not a member of the Cyg OB8 
association, but possibly a member of one of the other associations in that direction. 



Discussion 



The five systems studied here are all relatively long period, low semiamplitude SBls. The 
velocity semiamplitudes of these systems range from 13.1 km s _1 to 26.2 km s _1 with orbital 
periods ranging from 6.2 to 15.8 days. None of the se systems show eclipses or photometric 
variation in eit her the All Sky Aut omated Survey ( jPojmahskil 120021 ) nor The Amateur Sky 
Survey (TASS0; broege et al.ll2006t). T he only object with published time series photometry 



is HD 152147, andlSterken et al. 



(119931 ) f ound no var iability. BD— 16°4826 is an X-ray source 



but with a flux typical of single O-stars ( iNazel 120091 ) . 



4 http://sallman. tass-survey.org/servlet/markiv/ 
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For circular orbits, the orbital semiamplitude is given by 



Ki = 289 km s _1 (Mi/25M ) 1/3 (P/1O d)~ 1/3 - 



l + g)V3 



sin? 



where the mass ratio is q = MijM\. Therefore the low K\ systems in our sample are drawn 
from those with low mass companions (small q) and/or low inclination. W e can estimate 
the most probable q from the statistical method of iMazeh &: Goldberg (ll992f ). This method 
requires assumptions about the mass ratio distr i bution of the sample and the mass of the 
primary, which we estimate from iMartins et al.l (120051 ) . Based on observations of O-stars 
( jMason et al.lll998l ). we assume the mass ratio distribution is uniform. The resultant mass 
ratios (q) from this distribution are listed in Table HJ Estimates for the spectral classification 
of the secondary ar e made by comparing the mass in column 5 to the values listed in Table 
B.l of I Gray! ( 120051 ). We assume these companions are all main sequence stars, because all 
the primary stars are so young. Column 7 lists t he diff e rence in V magnitude between the 
primary and secondary, again as estimated from iGrayl (120051 ). These are all greater than 
2.3 mag, and they lead to the monochromatic flux ratios given in column 8. The extremely 
small flux ratios may be part of the reason why companions are not seen in the data. 

The true nature of the mass ratio distribution for O-type stars is unknown. (For a 
discussion of t he current state of our unders tanding; of the mass ratio distribution, see the 
recent work by lKiminki fc Kobulnickyl (120121 ) .) To explore the range of possibilities besides 
the flat distribution used above, we take two extreme values for a power law di stribution 
N(q) = N q a . The first is the standard Salpeter distribution ( ISalpeterl Il955l ) with an 
exponent of a = —1.35. The second has an exponent of a = 1.0 , a distribution that is 
more biased toward near equal- mass binaries, as was suggested in iPinsonneault fc Stanek 
(120061 ). With these mass ratio distributions, and estimates for the masses of the primary and 
secondary, we can solve for the probable inclination of the orbit based on the mass function. 
The resulting ranges of mass ratio and inclination are listed in Table |5j These ranges are 
not large, and the estimates for orbital inclination are low (i < 37°). 

Taking our value s for i and orbital period, combined with an estimate of the stellar radius 
( IMartins et al.ll2005l ). we compute the synchronous rotational velocity, v sinz (sync), listed 
in column 10 of Table HJ These numbers are all much smaller than the vsini measurements 
listed in column 5 of Table [31 which is consistent with the idea that these systems are too 
young to have attained synchronous rotation. 

Do the line blending effects of the companions influence our results for vsini or K{1 
To check the magnitude of the effect, we can consider the maximum separation of the spec- 
tral features of the individual components. The ratio of maximum Doppler separation to 
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projected rotational velocity is 

Ki(l + 1/g) 
f sin 2 

By using the mass ratio as computed above, we calculate this separation and list it in column 
9 of Table HI These values are all well below unity, meaning the lines of the components 
are always blended with those of the primary (except in cases where a lower ionization line 
appears in the secondary but not in the primary spectrum). The additional fact that the 
companions are extremely faint means that line blending problems related to measurements 
of v sin i and K\ will be minimal. As an example, we introduced a model spectrum cor- 
responding to a B3 V star, as estimated for the companion to BD — 16°4826 (Column 6, 
Table HJ), at a shift of K\{\ + 1/g) and normalized according to the flux ratio. This object 
was chosen because it has the largest ratio of Ki(l + l/q)/vsmi, so it should show the largest 
affect of blending problems. We chose a v sin % of the companion of 20 km s -1 . Such very nar- 
row lines would have more influence on the measured parameters of the primary than wider 
shallower lines. Measuring the velocity in the same way as described above, we compute a 
maximum offset of 3.6 km s _1 for velocity measurements, comparable to the measurement 
uncertainties on an individual spectrum. We similarly added a secondary model spectrum 
at ±Ki(l + 1/g) to each of the lines used to measure v sini, and calculated a value 5.0 km 
s _1 higher than previously. Therefore, even for the narrowest lined star in our sample, the 
uncertainty in vsini introduced by the companion is on the order of 5 km s -1 . Although the 
effects are small, the influence of a companion spectrum on the velocity measurements will 
be to pull the measurement back toward the rest wavelength of the line. This will make K\ 
lower, the mass function lower, and thus the estimate of mass ratio and orbital inclination 
lower. Consequently, the values of q and % in Table |5] are therefore probably lower than the 
true values. 

Companions will also have an influence on T e g and log g estimates. In the spectral 
range covered by our data, an O-star has several He II lines that a mid B-star does not. 
Therefore, the effect will be to make the O-star lines appear shallower, making the O-star 
appear cooler. He I lines are stronger in B-stars, and in combination with an O-star spectrum 
will make these lines appear deeper, again causing the O-star spectrum to appear cooler. 
Thus, our estimates for T e g and log g are also lower limits and are likely more indicative of 
the combined spectrum of the binary. 

Our results from this analysis suggest that all the objects in this sample are SB1. This 
is due to three main effects: low orbital inclination, broad spectral features of the primary 
stars, and faint companions. In the future, high resolution and high S / N spectra may b e able 



to extract the companion spectrum in systems like those studied here ( IGies et al.lll994f ). The 



next paper in this series will describe double-lined spectroscopic orbits for Galactic O-type 
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stars where companion spectra are clearly seen. 

We are grateful to the directors and staffs of KPNO and CTIO for their help in obtaining 
these observations. This material is based upon work supported by the National Science 
Foundation under Grant No. AST-0606861 and AST-1009080. This research has made use 
of the Washington Double Star Catalog maintained at the U.S. Naval Observatory. 

Facility: CTIO: 1.5m Facility: KPNO:2.1m 
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Table 1. Radial Velocity Measurements 



Star 


Date 


Orbital 


V r 


a 


O - 


C 


Name 


(HJD - 2,450,000) 


Phase 


(km s- 1 ) 


(km s- 1 ) 


(km f 




HDE 308813 


2714.805 


0.699 


-0.6 


1.7 




-2.2 


HDE 308813 


2715.708 


0.841 


20.6 


1.7 




4.6 


HDE 308813 


2715.840 


0.862 


31.2 


1.8 




8.2 


HDE 308813 


2716.707 


0.999 


27.3 


1.7 




-3.3 


HDE 308813 


2716.866 


0.024 


20.7 


1.9 




-9.7 


HDE 308813 


2717.697 


0.155 


32.8 


1.9 




12.0 


HDE 308813 


2717.801 


0.171 


16.3 


1.8 




-2.6 


HDE 308813 


2718.699 


0.313 


4.5 


1.8 




4.7 


HDE 308813 


3017.743 


0.484 


-13.6 


2.0 




0.0 


HDE 308813 


3018.753 


0.644 


-1.6 


1.8 




3.6 


HDE 308813 


3019.735 


0.799 


16.2 


2.0 




1.0 


HDE 308813 


3019.836 


0.815 


26.0 


1.9 




8.8 


HDE 308813 


3020.718 


0.954 


23.6 


1.9 




-6.1 


HDE 308813 


3021.713 


0.111 


28.1 


1.9 




2.6 


HDE 308813 


3021.816 


0.127 


29.6 


1.9 




5.7 


HDE 308813 


3022.721 


0.270 


-3.7 


1.9 




-9.5 


HDE 308813 


3022.865 


0.292 


-4.7 


1.9 




-7.4 


HDE 308813 


3023.702 


0.424 


-18.5 


2.0 




-7.3 


HDE 308813 


3023.806 


0.441 


-11.2 


1.9 




0.9 


HDE 308813 


3027.727 


0.059 


24.9 


2.0 




-4.2 


HDE 308813 


3028.704 


0.213 


16.4 


2.0 




3.0 


HDE 308813 


3028.787 


0.227 


19.6 


2.0 




7.8 


HDE 308813 


3029.744 


0.377 


-3.8 


2.0 




3.7 


HDE 308813 


3030.694 


0.527 


-11.2 


1.9 




2.1 


HDE 308813 


3030.776 


0.540 


-7.9 


2.0 




5.1 


HDE 308813 


3031.712 


0.688 


-1.4 


1.9 




-1.5 


HDE 308813 


3032.636 


0.834 


10.0 


1.9 




-9.6 


HDE 308813 


3032.782 


0.857 


24.5 


2.0 




2.2 


HDE 308813 


3033.713 


0.004 


34.9 


1.9 




4.3 


HDE 308813 


3159.469 


0.840 


15.0 


1.8 




-5.4 


HDE 308813 


3160.463 


0.997 


26.8 


1.8 




-3.8 


HDE 308813 


3162.531 


0.323 


-1.9 


1.8 




-0.5 


HDE 308813 


3163.532 


0.481 


-20.8 


1.9 




-7.3 


HD 152147 


3152.593 


0.668 


-46.7 


2.6 




-0.2 


HD 152147 


3153.620 


0.742 


-33.6 


2.5 




7.2 


HD 152147 


3154.580 


0.811 


-43.7 


2.5 




-8.6 


HD 152147 


3154.811 


0.828 


-35.1 


2.5 




-1.2 


HD 152147 


3155.743 


0.896 


-32.0 


2.7 




-2.3 


HD 152147 


3156.775 


0.970 


-23.8 


3.0 




3.5 


HD 152147 


3157.764 


0.042 


-26.7 


2.5 




0.7 


HD 152147 


3157.874 


0.050 


-40.4 


3.4 




-12.7 


HD 152147 


3158.696 


0.109 


-29.2 


2.5 




0.8 


HD 152147 


3159.714 


0.183 


-35.0 


2.4 




-0.2 


HD 152147 


3160.721 


0.256 


-41.6 


2.5 




-1.0 


HD 152147 


3161.737 


0.329 


-39.2 


2.5 




7.1 
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Table 1 — Continued 



Star 


Date 


Orbital 


Vr 


a 


O — C 


Name 


(HJD — 2,450,000] 


Phase 


(km s _1 ) 


(km s _1 ) 


(km s~ x ) 


HD 152147 


3162.688 


0.398 


-60.2 


2.6 


-9.5 


HD 152147 


3163.709 


0.472 


-50.8 


2.6 


2.2 


HD 152147 


3164.599 


0.536 


-57.7 


2.7 


-4.8 


HD 164536 


31 52 7^3 


0.559 


—26.6 


2.0 


2.7 


HD 164536 


31 53 777 


0.636 


-27.4 


1.9 


-5.3 


HD 164536 


3153 916 


0.646 


-17.0 


1.9 


3.7 


HD 164536 


3154 781 


0.711 


-7.5 


1.8 


3.8 


HD 164536 


3155.720 


0.782 


2.1 


1.9 


1.9 


HD 164536 


31 56 633 


0.850 


14.4 


2.0 


3.9 


HD 164536 


3156 910 


0.871 


4.8 


2.0 


-8.3 


HD 164536 


3157.881 


0.943 


15.1 


2.0 


-4.5 


HD 164536 


3158.803 


0.012 


20.4 


1.9 


-0.7 


HD 164536 


3159.750 


0.083 


21.0 


1.9 


3.2 


HD 164536 


3160.773 


0.160 


19.0 


1.9 


9.9 


HD 164536 


3161.842 


0.240 


-4.9 


2.1 


-1.5 


HD 164536 


3162.789 


0.311 


-26.5 


2.0 


-11.8 


HD 164536 


3163.779 


0.385 


-19.5 


1.9 


5.0 


HD 164536 


3164.751 


0.458 


-31.3 


2.0 


-1.1 


HD 164536 


3164.870 


0.467 


-30.8 


2.0 


-0.3 


BD-16°4826 


3152.899 


0.391 


-0.7 


2.5 


-0.5 


BD-16°4826 


3153.797 


0.445 


2.9 


2.3 


5.4 


BD-16°4826 


3153.929 


0.453 


1.6 


2.5 


4.3 


BD-16°4826 


3154.913 


0.513 


-6.1 


2.3 


-2.8 


BD-16°4826 


3156.646 


0.617 


3.9 


2.4 


3.5 


BD-16°4826 


3157.914 


0.693 


5.2 


2.4 


-0.7 


BD-16°4826 


3158.901 


0.753 


13.8 


2.3 


2.6 


BD-16°4826 


3159.902 


0.813 


13.4 


2.5 


-2.9 


BD-16°4826 


3161.900 


0.934 


26.6 


2.7 


2.9 


BD-16°4826 


3162.901 


0.994 


22.3 


2.4 


-2.6 


BD-16°4826 


3164.888 


0.114 


20.2 


2.4 


-1.2 


HDE 229232 


4785.649 


0.198 


-37.5 


2.5 


4.0 


HDE 229232 


4786.594 


0.351 


-56.1 


2.2 


-0.2 


HDE 229232 


4787.576 


0.511 


-69.7 


2.1 


-7.5 


HDE 229232 


4788.576 


0.673 


-47.4 


2.2 


6.5 


HDE 229232 


4789.583 


0.836 


-37.1 


2.1 


1.4 


HDE 229232 


4790.572 


0.997 


-37.5 


2.1 


-6.5 


HDE 229232 


4791.574 


0.160 


-36.9 


2.1 


1.3 
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Table 2. Circular Orbital Elements 



Parameter 


HDE 308813 


HD 152147 


HD 164536 


BD-16°4826 


HDE 229232 


P(d) 


6.340±0.004 


13.8194±0.0003 


13.4±0.6 


15.8±1.3 


6.2±0.2 


T (HJD - 2,400,000.0) 


53033.69±0.08 


31212.0±0.2 


53158.6±0.8 


53162. 8±1. 8 


54790.6±0.1 


7 (km s _1 ) 


8.5±1.1 


-40.1±0.7 


-4.9±1.1 


11. Oil. 


-46.6±0.8 


Ki (km s" 1 ) 


22.2±1.5 


13.1±1.0 


26.2±1.5 


13.1±1.6 


15.6±1.2 


/(m) (M ) 


0.007±0.001 


0.0032±0.0007 


0.025±0.005 


0.004±0.002 


0.002±0.001 


ai sin i (-Rq) 


2.8±0.2 


3.6±0.3 


6.9±0.5 


4.1±0.6 


1.9±0.2 


rms (km s _1 ) 


6.1 


5.7 


6.1 


3.2 


7.3 


N 


33 


31 


16 


11 


7 



Table 3. Individual Stellar Parameters 





Spectral 




log a 


v sin i 


Spectral Classification 


M v 


E(B - V) 


Distance 


Object Name 


Classification 


(kK) 


(cm s~ 2 ) 


(km s _1 ) 


from (T eff ,log g) 


(mag) 


(mag) 


(kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HDE 308813 


09.5 V 


29.9±0.3 


3.73±0.09 


204±10 


B0 V a 


-3.90 


0.30 


2.82 


HD 152147 


09.7 lb 


27.8±0.5 


3.10±0.06 


150±28 


09.5 I a 


-6.28 


0.61 


2.11 


HD 164536 


09 III 


37.4±0.9 


4.25±0.17 


230±14 


06.5 V b 


-4.77 


0.22 


1.78 


BD-16°4826 


05 


39.9±6.3 


4.04±0.40 


131±28 


05.5 V b 


-5.07 


1.24 


1.60 


HDE 229232 


04 Vn((f)) 


41.7±1.3 


4.05±0.14 


273±19 


05 V b 


-5.21 


1.12 


1.79 



Calibration of iBohm-Vitense! Jl98ll) . 
b Calibration of lMartins et al.1 d2005h . 



Note. — Spectral classifi cation referen ces: HDE 308813 : ISchildl lll970h ; HP 152147: | Walbornl lll972h : HD 164536: 
iMacConnell fc Bidelmanl 1119761) ; BD-16°4826: iHiltner fc Johnson! lll956h ; HDE 229232: IWalbornl lll97aV 
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Table 4. Probable Unseen Companion Parameters 





f(m) 


Mi 




M 2 


Estimated Secondary 


Ay 


Flux Ratio 


K^l + l/q) 


vs'mi (sync) 


System 


(M ) 


(Mq) 


1 




Spectral Classification 


(mag) 


(/2//1) 


vs'mi 


(km s- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


HDE 308813 


0.007 


14 


0.26 


3.4 


B7 V 


2.7 


0.08 


0.53 


20 


HD 152147 


0.003 


28 


0.18 


5.0 


B4 V 


4.8 


0.01 


0.55 


23 


HD 164536 


0.025 


28 


0.30 


8.3 


B2 V 


2.3 


0.12 


0.49 


13 


BD-16°4826 


0.004 


34 


0.19 


6.4 


B3 V 


3.1 


0.06 


0.67 


9 


HDE 229232 


0.002 


38 


0.16 


5.9 


B3 V 


3.2 


0.05 


0.45 


29 



Table 5. Mass Ratios and Inclinations for Various Mass Ratio Distributions 







q 










System 


a = -1.35 


a = 0.0 


a = 1.0 


a = -1.35 


a = 0.0 


a = 1.0 


HDE 308813 


0.15 


0.26 


0.43 


35.5 


20.5 


6.0 


HD 152147 


0.10 


0.19 


0.38 


30.4 


16.6 


10.5 


HD 164536 


0.18 


0.30 


0.45 


37.0 


22.7 


12.6 


BD-16°4826 


0.11 


0.19 


0.39 


32.4 


18.3 


13.1 


HDE 229232 


0.09 


0.16 


0.36 


28.0 


15.4 


7.4 
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Fig. 1. — Radial velocities for HDE 308813. Filled circles are data f rom this analysis wit h 
uncertainties shown by line segments, and open circles are data from lHuang fc Giesl ( 120061 ). 
Plotted beneath the velocities are the observed minus calculated values based on the circular 
orbital fit. Also shown is the 5 km s _1 uncertainty associated with the wavelength calibration. 
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Fig. 2. — Radial v elocities for HP 15214 7. Filled circles are dat a from t his an alysis, the 
open square is fro mlConti fc Ebbets open triangles are from IStruvd ( 1l944l ). and open 

f 19881 ). with uncertainties shown by line segments. Plotted 



diamonds are from Levato et al. 



beneath the velocities are the observed minus calculated values based on the circular orbital 
fit. 
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Fig. 3. — Radial velocities for HD 164536. Filled circles are data from this analysis, with 
uncertainties shown as line segments. Plotted beneath the velocities are the observed minus 
calculated values based on the circular orbital fit. 
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Fig. 4. — Radial velocities for BD— 16°4826. Filled circles are data analyzed here, with 
uncertainties shown as line segments. Plotted beneath the velocities are the observed minus 
calculated values based on the circular orbital fit. 
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Fig. 5. — Radial velocities for HDE 229232. Filled circles are data from this analysis, with 
uncertainties shown as line segments. Plotted beneath the velocities are the observed minus 
calculated values based on the circular orbital fit. 



-26 - 









o 

CO O 

3 x - 



X 



5 s 

3 -t 



CD 

oo 

CD 



ca 
x 



CD 
X 



CD 
X 



CD CD CD 
XX X 

J I I I I I L_ 



-i HDE 229232 
(05 V) 

BD-16°4826 
(05.5 V) 

HD 164536 
(06.5 V) 

HDE 308813 
(BO V) 

HD 152147 
(09.5 I) 



4000 4200 



4400 4600 4800 5000 
Wavelength (Angstroms) 



5200 5400 



Fig. 6. — Spectra and prominent spectral features for objects discussed here. The hottest 
star is at top and targets are then plotted in order of decreasing effective temperature. 



